Spin Dynamics of NiCl 2 -4SC(NH 2 ) 2 in the Field-Induced Ordered Phase 
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NiCl2-4SC(NH2)2 (known as DTN) is a spin-1 material with a strong single-ion anisotropy that is regarded 
as a new candidate for Bose-Einstein condensation (BEC) of spin degrees of freedom. We present a systematic 
study of the low-energy excitation spectrum of DTN in the field-induced magnetically ordered phase by means 
of high-field electron spin resonance measurements at temperatures down to 0.45 K. We argue that two gapped 
modes observed in the experiment can be consistently interpreted within a four-sublattice antiferromagnet model 
with a finite interaction between two tetragonal subsystems and unbroken axial symmetry. The latter is crucial 
for the interpretation of the field-induced ordering in DTN in terms of BEC. 

PACS numbers: 75.40.Gb, 76.30.-v, 75.10.Jm 
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Introduction.- Field-induced phase transitions in magnets 
have recently received a considerable amount of attention 
id H S H S S 01, particularly in the context of the so- 
called Bose-Einstein condensation (BEC) of spin degrees of 
freedom. In accordance to the BEC scenario, for a uniform 
gas of identical particles, at some finite temperature T c when 
the de Broglie wavelength becomes comparable to the aver- 
age distance between the particles, a macroscopic fraction of 
the gas can be "condensed" into a single coherent quantum 
state and BEC occurs. It is worth mentioning that although 
there are some important arguments |@1 about the interpre- 
tation of the field-induced transitions in magnets in terms of 
BEC (based on its original, canonical definition |9J]), the ap- 
plication of this formalism to grand-canonical ensembles of 
quasi-particles in magnets appears to be widely accepted. An 
important property of BEC is the presence of U(l) symmetry, 
which corresponds to the global rotational symmetry of the 
bosonic field phase. Below T c , the U(l) symmetry sponta- 
neously gets broken, the wave function of the condensate gets 
coherent on a macroscopic scale and, as a consequence, a gap- 
less Goldstone mode is acquired. The model of BEC in mag- 
nets assumes that the spin-Hamiltonian is axially symmetric 
with respect to the magnetic field, requiring U(l) rotational 
symmetry above T c . 

In accordance with mean-field BEC theory, the phase- 
diagram boundary for a three-dimensional magnet should 
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obey a power-law dependence, H — H c i ~ T c 1 , where 
H c i is the first critical field as T — > 0. Importantly, criti- 
cal exponents extracted from phase diagrams alone can not 
be regarded as sufficient criteria for identifying field-induced 
transitions as BEC. For instance, the re-opening of the energy 
gap in the excitation spectrum of TlCuCi3 (which, based on 
analysis of critical exponents, was regarded as the best real- 
ization of BEC of spin degrees of freedom in magnets J^) in 
the field-induced ordered state [10] is a clear evidence for a 
broken uniaxial symmetry, which rules out the description of 
the magnetic ordering in this compound in terms of BEC. 



NiCl 2 -4SC(NH 2 ) 2 (known as DTN) is a gapped S = 1 
system with a single-ion anisotropy D dominating over the 
exchange coupling J Eta , and is a new candidate for BEC 
of spin degrees of freedom. Although below T c < 1.2 K 
DTN exhibits a field-induced antiferromagnetic (AF) order- 
ing, with critical fields B c \ = 2.1 T and B C 2 = 12.6 T 
and exponents corresponding to the BEC scenario 10], the 
low-energy excitation spectrum in DTN (and, correspond- 
ingly, the microscopical picture of magnetic interactions) in 
the field-induced ordered state still remains an open question. 
In the present work, we have addressed these issues exper- 
imentally by means of electron spin resonance (ESR) mea- 
surements performed at temperatures down to 0.45 K. Two 
gapped modes were observed. Based on our detailed analysis, 
it is argued that the above excitation spectrum can be consis- 
tently interpreted within a four-sublattice AF model with an 
intact axial symmetry (at least on the energy scale down to 
1.2 K, which corresponds to the lowest frequency used in our 
experiments, 25 GHz). The latter is of particular importance, 
being a necessary prerequisite for the interpreting of the AF 
ordering in DTN in terms of the BEC scenario. 

DTN is characterized by the 14 space group lfl2ll with a 
body-centered tetragonal lattice that may be viewed as two 
interpenetrating tetragonal subsystems (hereafter TS). At B \\ 
c the spin dynamics can be described by the spin-Hamiltonian 
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where S n are spin- 1 operators at site n, the vectors d connect 
the site n to its nearest neighbors within the same subsys- 
tem, h = g c HBB is the Zeeman term, and 7Yi n t describes (yet 
unspecified) additional interactions. Assuming that the latter 
are much weaker than the interchain interaction within TS, 
the Hamiltonian parameters were estimated as D = 8.9 K, 
J c = 2.2 K, J a . b = 0.18 K, and g c = 2.26 |l3]. 
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FIG. 1: Typical ESR spectra of DTN in the quantum-disordered 
phase taken at T = 1.4 K (left) and in the ordered phase taken at 
T = 0.45 K (right). The spectra are offset and rescaled for clarity. 
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FIG. 2: Temperature dependencies of ESR spectra taken at v = 
143.2 GHz (left) and 28.8 GHz (right); the absorption scale is re- 
tained at each frequency. The inset shows the temperature evolution 
of the energy gap Aa" of the K mode, at magnetic field of 8 T (the 
dashed line is a guide to the eye). 



Experimental results- The ESR measurements were 
done at the Kapitza Institute using a transmission-type ESR 
spectrometer equipped with a cylindrical multimode resonator 
and a 3 He-cryostat. High-quality single-crystalline DTN sam- 
ples from the same batch as in Ref. 1 13] were used. The mag- 
netic field was applied along the tetragonal c axis. 

The examples of the ESR spectra taken at 1.4 K are shown 
in Fig. Q] (left), while the ESR temperature evolution is illus- 
trated in Fig. |2] The corresponding frequency-field depen- 
dencies of magnetic excitations are presented in Fig. [5] Two 
ESR lines (denoted as A and C) have been observed at 1.4 K. 
Comparison with the data of Ref. 1 1311 (filled grey circles, Fig. 
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FIG 3: (Color online) The frequency-field dependence of ESR in 
DTN measured at T = 1.4 K (open circles) and 0.45 K (squares and 
stars). Filled grey circles denote the high frequency data (T — 1.6 K) 
taken from Ref. fiUl . The dashed (red) lines correspond to results 
of calculations for the simplest axially symmetric two-sublattice AF 
model with parameters g c = 2.26, D = 9.4 K, J = 2.0 K. The 
solid (cyan) lines correspond to results of model calculations assum- 
ing d = 0.02 K (see text for details). 



reveals that the mode A continues smoothly from the ESR 
line in the low-field disordered phase, while the mode C cor- 
responds to single-magnon excitations in the high-field phase. 

Several important features have been observed at lower 
temperatures. Upon cooling down, the modes A and C 
shift towards each other (Fig. O. While at a temperature of 
1.4 K they seem to cross at zero frequency in the vicinity of 
B ~ 8 T, in the low-temperature AF ordered phase they are 
converted into a new gapped mode K. This mode exhibits 
a slight but distinct splitting, the origin of which will be dis- 
cussed below (the corresponding frequency-field dependence 
of the resonances is denoted in Fig. [3] by pairs of open and 
closed squares). The temperature dependence of the gap Ak 
at B = 8 T has been obtained by measuring the temperature of 
the maximum absorption at a fixed field (i.e., by the observa- 
tion of so-called "temperature resonance") or by using a con- 
ventional ESR procedure, recording transmission vs magnetic 
field at different temperatures (circles and squares in Fig. [2] 
inset, correspondingly). The diminishing of Ak upon warm- 
ing can be explained by the decrease of the order parameter. 
The low temperature spectrum demonstrates one more set of 
resonances appearing at low frequencies (mode L, denoted by 
stars in Fig. [3]). Unlike mode K, the integrated intensity of 
line L is 50 — 100 times smaller than that of modes outside 
the AF phase, e.g. modes A and C. 

Discussion- The observation of the gapped mode L and 
splitting of mode K in a field range of 6 -j- 10 T clearly indi- 
cates the presence of additional interactions, that are not ac- 
counted for in the simplest axially symmetric two-sublattice 
AF model (which would have the Goldstone mode as its low- 
est energy excitation and a single mode K as the next excita- 
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tion branch as indicated by the dashed (red) lines in Fig . [3J ■ To 
gain further insight, we have studied several possible mecha- 
nisms that might cause gaps in the ESR excitation spectrum. 

To do so, we have considered the model (JTJ with various 
additional interactions Hint, using the mean-field approach 
outlined below. Assuming the wave function to be a prod- 
uct of single-spin coherent states parameterized by a pair of 
two-component real vectors u, v: 



|V>) = s\to) + ^ ( U 3 + ^i)!^'}: 
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where the spin-1 states are |0) = \to) and |±) = =p-^=(|t x ) ± 
i\t y )), one obtains the effective Lagrangian in the form C = 

h J2n( U n ■ dtV n - V n ■ d t U n ) - (U) , with 



(W> = E n {D(ui + vi)-2h-(u n xv n )} 

+ 2 J2(nn>) Jn-n>{(Vn ■ V n > ) (s n S n > + U 

){v 

int/ i 
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where (nn 1 ) denotes summation over nearest neighbors 
within the corresponding TS. For weak fields, the static mean- 
field solution is u, v = 0. However, at B > B c \ one obtains a 
nontrivial solution with finite u and v determined by the min- 
imization of (TL), which describes the in-plane AF order and 
finite magnetization along the z axis. Linearizing around the 
static solution and solving the equations of motion, one can 
find the excitation energies Outside the AF phase the ESR 
modes correspond to Ek=o, while inside the ordered phase, 
due to the folding of the Brillouin zone, one has the additional 
transitions at Qb = (ft/a, n/b, 7r/c) with energies £k=Q B - 

First we have considered the possibility of axial symmetry 
breaking, caused either by rhombic in-plane anisotropy or by 
the Dzyaloshinskii-Moriya (DM) interaction inside a TS, with 
the DM vector d deviating from the anisotropy c axis (note 
that DM is generally allowed in DTN due to the absence of 
inversion symmetry [14]). In both cases, the axial symmetry 
break-down would open a gap in the excitation spectrum, 'lift- 
ing' the Goldstone mode. Even so, this mode corresponds to 
a coherent collective excitation of the magnon condensate be- 
low T c , and thus would be expected to be particularly strong in 
comparison with ESR excitations in the disordered phase. As 
mentioned, the gapped mode L observed in our experiments 
is much less intense than ESR absorptions above T c , raising 
questions about the above scenarios. In addition, the rhombic 
in-plane anisotropy would necessarily cause splitting of the 
doublet at B = (roughly of the same strength as the max- 
imum energy of the L mode), which was not observed in the 
experiment. It is worth mentioning that the DM term is cubic 
in u, v, so it does not contribute to the quadratic spectrum at 
B < B c i and thus does not produce mode splitting at B = 0. 
However, a finite in-plane component of d (e.g., d y ) con- 
tributes to the energy as Wdm = \ Yl n s r lnd y {S n x S n +s) y , 
where r\ n = ±1 for sites n belonging to the different mag- 
netic sublattices (inside the same TS). Close to the fully spin- 
polarized state, this contribution provides an energy gain lin- 



ear in the AF order parameter, while the corresponding Zee- 
man energy loss is quadratic, so AF order is favored for arbi- 
trarily large B. Consequently (and this is of particular impor- 
tance for our analysis), such a symmetry-breaking DM term 
would lead to the absence of the second critical field B C 2, 
which is incompatible with results of flO]. 

Here we argue that a natural explanation of all available ex- 
perimental results is possible if we take into account a weak 
isotropic corner-center interaction of magnetic ions in the 
body-centered tetragonal lattice, i.e. an interaction between 
TS that preserves the axial symmetry. In case of such interac- 
tions in addition to conventional (relativistic) modes, modes 
with antiphase oscillations of two interacting AF sublattices 
(exchange modes) should be present. Having much weaker 
coupling to the microwave field 11511 . exchange modes should 
be less intensive than relativistic modes. In accordance to our 
observations, the gapped mode L would then be an exchange 
mode, coexisting with the lowest-energy relativistic mode. 

The physically simplest scenario would correspond to the 
isotropic "corner-center" exchange. However, the theoretical 
analysis of such a model is very difficult; if, as indicated by 
previous studies [13], the exchange interaction within each TS 
is AF, then the system is highly frustrated and its mean-field 
ground state at B > B c \ is infinitely degenerate. This degen- 
eracy may be lifted by quantum fluctuations or by additional 
interactions. For that reason, we would like to illustrate the 
effect of corner-center coupling by assuming a finite DM in- 
teraction between the TS, which already lifts the degeneracy 
at the mean-field level: 



Hint = ^ Vnlln'd • (S A ,n X S 



B,n' 



(3) 



where the DM vector d is along the c axis and A, B denote 
the two TS. Such interaction favors a 90° angle between the 
AF order parameters of the subsystems. The static mean-field 
solution at B > B c \ is then 



VA,t 



VB,r 



= Vr. 



u x u y 
v y v x 
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and the excitation energies are obtained from the secular 

r f u 

equation det(i?) = 0, with R = TT+ ,_, , where 
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the matrix U describes intersubsystem interaction, 



U = Ad,, 





si - 2v 2 - 
-u v (l - 




vl/sl) 





-u v (l - vl/sl) 
ulvl/sl 




4 



and the following shorthand notations have been used: 

Wi = -4s§J fc + w 9 , W 2 = u v Q (f k - g k ) 

W 3 = 4J k (2u 2 - 1) + i>g(3/ fc - g k ), W A = 4u Q v Q Jo 

W 5 = +Vq f k - u%g h , W 6 = -4^(J + d ) 

with Jfc = J2\= a bc^e cos(fc • A), do = 4d, and 
d k = (d /8) JJ (l — exp(ife • A)) 

A— a,b,c 

/fe = 4(J fc - J - do), 9k = {4vl/sl)(J k + Jq + d ). 

The results of model calculations using g c = 2.26, D — 
9.4 K, J = 2.0 K, d = 0.02 K are shown in Fig. |3]by solid 
(cyan) lines. Due to a finite interaction between the two TS 
(the term d), the low-energy mode is split into a doublet with 
a gapped upper component and a zero frequency lower com- 
ponent. One can see that the model qualitatively describes 
frequency-field dependencies of all observed ESR modes as- 
suming the existence of a gapless Goldstone mode G, which 
can not be detected experimentally. It is important to mention 
that it was not possible to fit the frequency-field dependence of 
magnetic excitations in DTN in the AF phase using the set of 
parameters obtained in Ref. 11311 . although they are very close 
to those used in our calculations. This discrepancy mainly 
stems from neglecting quantum fluctuations in the mean-field- 
theory approach used in the present paper. In the above cal- 
culation the mode K (determined by £fc = o) comes out doubly 
degenerate, which is an artefact of the model assumption on 
purely DM inter-TS interaction, making the interaction matrix 
U vanish at k — 0. Any finite symmetric exchange interaction 
between the TS will be sufficient to lift this degeneracy and 
thus will explain the observed slight splitting of the K mode. 
The theory also predict the existence of a third ESR mode 
M in the AF phase (Fig.|3), which corresponds to the second 
magnon branch at fe = Qb ', this mode could not be observed 
in the present study due to the limited frequency range. 

In summary, high-field ESR studies of magnetic excitations 
in the field-induced ordered phase of DTN have been per- 



formed at frequencies down to 25 GHz. Two gapped modes 
were observed in the ESR spectrum. Our experimental ob- 
servations can be consistently interpreted within the four- 
sublattice AF model with an intact axial symmetry, which is of 
crucial importance for the interpretation of the field-induced 
ordering in DTN in terms of BEC of spin degrees of freedom. 
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